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Abstract 

Motivated by ideas from the consistent histories approach to quan- 
tum mechanics, we examine a simple model of hydrodynamic coarse 
graining for a scalar field. It consists in averaging the field over spa- 
tial regions of size L and constructing the evolution equations for the 
coarse grained quantities, identifying dissipation and noise. 
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1 Introduction 



1.1 Motivation 

In modern cosmology there seems to exist a gap between our description 
of early universe using essentially quantum field theory (hence potentially 
involving highly non-classical behaviour of matter) and the later time (up 
to current era) description of matter as essentially classical by a few phe- 
nomenological parameters (as corresponding for instance to a classical fluid). 

This gap is not conceptual: one can intuitively visualise that as the scale 
factor becomes larger somehow the quantum behaviour of the fields gets 
smoothed out and matter falling to lower energies can essentially be de- 
scribed by a collective set of variables, (usually hydro dynamic). The prob- 
lem arises when one tries to give a more precise description of the underlying 
mechanisms of this process. 

It can be safely said that the quantum to classical transition for quan- 
tum field theories has not been understood very well yet. This is in spite 
of the large amount of activity in recent years and the application of rela- 
tively recent concepts and ideas, successful in the realm of non relativistic 
quantum mechanics (the environment superselection rules approach , con- 
sistent /decoherent histories [0, |, ||). 

The issue is more strongly highlighted in the context of inflationary mod- 
els; a necessary ingredient in them is the amplification of quantum fluctua- 
tions to classical so that the inhomogeneity can be properly explained. But 
so far there is no consensus on how this classicalisation is taking place and 
more importantly why. 

The problem of course stems from the extremely complicated nature of 
quantum field theories. Being a systems with an infinite number of degrees 
of freedom there is a large (both technical and conceptual ) difficulty into 
implementing schemes succesful in ordinary quantum mechanics. 

What is usually considered in such discussions is the phase space classi- 
cal limit of a quantum field, that is the emergence of a classical field limit 
from the full quantum field theory. In such a limit the expectation value 
of the field evolves according to the Hamilton equations plus the action of 
some stochastic terms. Since in such discussions usually free fields are con- 
sidered, the discussion essentially comes into the consideration of each mode 
separately. 
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The question of course exists whether the results coming from mode by 
mode sphtting are robust enough when one takes into account field inter- 
action. But even, letting this issue unanswered, there are is a number of 
problems. It is sometimes stated that modes become classical solely by the 
squeezing induced throught the time variation of the scale factor. There are 
strong arguments against this interpretation: from a phase space perspective 
sneezing entails a highly non - classical behaviour (see for a more thorough 
discussion) 

Inspired by the success of the environment induced decoherence program 
the coupling of the field to some kind of environment has been proposed as 
a mechanism for classicality. There are some problems with this approach 
two: is there a natural splitting between system and environment in early 
universe? decoherence in non-relativistic quantum mechanics appears only 
when environment is much larger than the system; is this the case in field 
theory or spaecial initial conditions are to be assumed? Another idea pro- 
fessed in this context is that the role of the environment could be played by 
the higher order (inaccesible to the local observer ) n - point funtion sector 
of the field 

This is an interesting possibility, which forms one of the two main mo- 
tivating ideas of this paper. The other (and main one) is the Gell - Mann 
and Hartle decoherent histories program [Q. In this formalism, the classical 
limit is obtained by the study of coarse grained temporal propositions (his- 
tories) about the physical system. Classicality is equivalent to a consistency 
condition between pairs of histories, quantified by the vanishing of the off - 
diagonal terms of the decoherence functional. 

Concerning our aims, the great strength and versatility of this approach 
lies in the fact that it provides a generic algorithm to determine whether a set 
of histories effectively becomes classical or not, irrespective of the physical 
quantities to which they refer. Hence, in principle, it is sufficient to specify 
the suitable coarse grained operators, to be able to determine the existence 
of a corresponding quasiclassical domain. The discussion of classicality is 
therefore not restricted to the phase space properties of the system. 

As we argued, in the case of field theory phase space clasicality is a 
rather difficult issue to study. But in order to obtain an effective, cosmologi- 
cally relevant description it would be sufficient to concentrate on much more 
coarse - grained observables: essentially hydrodynamic or thermodynamic 
ones. Indeed, the emergence of hydrodynamic behaviour in many body sys- 
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terns (starting from first principles and no assumptions as local equilibrium) 
has been persistently advocated by Gell-Mann and Hartle. The issue is very 
important (not only in the cosmological setting), but involves such a degree 
of technical complexity, that only very simple systems have been examined 
so far 0. 

In any case , we believe that concentration on hydrodynamic properties of 
quantum fields is an important and not sufficiently addressed avenue towards 
understanding the emergence of classicality in quantum field theory. A great 
potential advantage of this approach, is that the observables involved are 
much more coarse-grained than the phase space ones and hence classicality 
might be emergent even if the field does not have a good classical field limit. 

This paper We are not in a position to give a full treatment of the 
emergence of hydrodynamic behaviour in field theory, within or without the 
consistent histories framework. Our aim in this paper is much more modest: 
in the simplest possible field theory (scalar field in Minkowski spacetime) 
we shall try to see how a simple hydrodynamic type of coarse graining can 
be implemented, focusing mainly on the general quantitative features of the 
evolution of coarse grained quantities. We should stress here, that even in 
such a simple system the behaviour is anything but trivial. 

But before explaining our model in more detail, let us first give some 
general considerations. 

Any hydrodynamic description of a theory, whose full dynamics are known 
involves three spatial scales, a clean separation of which is usually necessary 
to render the description be succesful . The first one, Imic is the time-scale 
within which microscopic processes are taking place. This is essentially deter- 
mined both by the dynamics and the initial state of the system. The second 
scale, lav is the scale, within which microscopic processes are averaged out. 
A volume of the size of can be thought as being sufficiently and indepen- 
dently described by a spatial average of all relevant physical quanties. This 
scale is determined primarily by the previous history of the system, and at 
least in the field theory case its existence is thought to necessitate a particu- 
lar class of initial states. It is clear that we should have Imic << lav Finally, 
the third scale lobs corresponds to the level of observation. It is fixed by 
the external constraints to the system . If one wants to take the continuum 
approximation, one necessarily has to demand that lav « lobs, but even if 
this is not true an effective lattice type of description can be implementd for 
the coarse grained dynamics. 
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These consideations fix more or less the strategy we shall adopt in this 
paper. We are going to consider a coarse-graining operation for the field 
consisting in spatial averaging within a spatial region of volume L^. That 
is, we consider a foliation of spacetime and to each Cauchy hypersurface we 
construct a cubic lattice of side L (which is to be identified with lav) ■ The 
field averaged values are to be thought as our coarse grained variables and an 
effective equation giving the evolution of their expectation values is going to 
be constructed through the use of the Zwanzig projector technique P, |ll|, |1^]. 
This is in spirit similar to a technique frequently employed in astrophysical 
discussions. 

A number of remarks are necessary at this point. First of all our choice 
for coarse- graining is not Lorenz invariant. Clearly, a different choice of fo- 
liation would lead to different value for L and hence to different expressions 
for the evolution equations. In particular the condition Imic « lav niight 
not hold in a new Lorentz frame. It is in this sense that both our choice of 
foliation and the value of L is not to be thought asarbitrary, but rather as 
phenomenological, and contingent upon a particular (cosmological) history. 
This is to be understood in the following way: Suppose that at some stage 
of the cosmological evolution, we observe that there is a preferred class of 
hypersurfaces (possibly the surfaces of approximate isotropy), such that: 

1. Histories of field variables smeared in regions of volume decohere 

2. Imic as can be read from the dynamics and the initial state (essentially 
would correspond to the maximum of characteristic length scales in the de- 
coherence functional) is much smaller than L. 

3. Our scale of observation (which is to be thought in the order of magnitude 
of the scale factor) is much larger than L. 

Then the evolution equations for this coarse-grained field variables, can es- 
sentially be constructed along the lines we shall develop in this paper. We 
cannot help overemphasizing, that the value for L is not arbitrary. Rather 
we should consider it to correspond to the minimum value for which histo- 
ries of averaged field variables decohere, and should be computed from first 
principles did we have a satisfactory (and computationally tractable) history 
version of field theory (or any other classicalisation scheme) as well as a good 
specification of the initial condition 

We believe, that an analogy with low energy physics wold help to clarify 
the above discussion. Consider a box ( whose side we take to correspond to 
lobs)i in which at one paricular moment of time, we insert a large number of 
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nuclei and electrons. The length scale of interaction for this highly active gas, 
should be decided by weighting a mean of the de-Broglie wavelengths of the 
particles ( which is determined by the initial state) with some possible length- 
scale determined by electromagnetic screening (hence the Hamiltonian is also 
contributing). Given sufficient time (and assuming low initial energy), our 
experience says that clusters will form (atoms) the evolution of which will be 
substantially less "violent" . Now it may happen as we focus our attention to 
a particular region of the box (with volume L^), the average value of physical 
quantities, will change smoothly over time and will not exhibit irregularities 
of quantum mechanical nature. Hence to these averaged variables we shall 
be able to assign probabilities. Taking the smallest value for L for which this 
is possible, we shall have the possibility of giving a description of the system 
in terms only of such coarse-grained observables. If in addition, L turns out 
to be much smaller than the size of the box, then a continuity assumption 
can be invoked and we would finally arrive at a hydrodynamic description of 
the system. For this regime, in principle, we should be able to construct a 
Navier-Stokes type of equation. 

As we said we shall restrict ourselves to a free field case. Independently 
of the above considerations, this is of interest on its own right: in recent 
work on stochastic infiation [|10| coarse graining of the infiaton field has been 
considered corresponding essentially to smearing the field in horizon scale 
spatial volume. In such coarse graining, it was assumed that only the low 
energy modes significantly contribute. In our model we shall be able to verify 
the truth of such an assertion, as well identify a source of noise due to the long 
range correlations of the initial state which is overlooked in such discussions. 

We should add that actual decoherence of field variables seems to ne- 
cessitate interaction terms, but as soon as this is established, the quadratic 
part of the Hamiltonian is the one giving the dominant contribution to the 
hydrodynamic equations (at least in the perturbation theory regime). The 
reason for this is that the coarse graining operation does not commute with 
the Hamiltonian, and even if interaction terms are assumed, their contribu- 
tion is going to be of the next order in perturbation theory. In any case, 
there is no conceptual or technical difficulty in considering interaction terms 
( this can be done by the conjuction of our coarse graining operator with the 
one of §). 
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2 Implementing the coarse - graining 



2.1 The formalism 

The Zwanzig method As we mentioned in the introduction , we are going 
to employ the Zwanzig projection formahsm for the implementation of the 
coarse - graining |n|, |T^. We remind the reader that this consists in 
introducing an indempotent "coarse graining operator" in the space of states 

p^Prel = Pp P' = P (2.1) 

which projects the state in the level of description. The irrelevant part of 
the state is then given by 

Pirr = (1 - P)P (2.2) 

By duality through the trace functional a corresponding operator can be 
introduced in the space of observables and we can obtain the evolution equa- 
tion for the expectation values of relevant observbles i.e. the ones for which 
PA = y4 ( P is assumed self - adjoint) 

^^^{A){t) - (PLA){t) 

+i f dr (PL{1-P)e'^^{l~ P)LP A) {t-r) = F^it) (2.3) 
Jo 

where hA = [H, A] in terms of the Hamiltonian H and e~*^* is the corre- 
sponding evolution operator. 

FA{t) is a stochastic in nature term, since it depends on the irrelevant 
components of the initial state that are inaccesible from our level of descrip- 
tion. It reads 

F^(t) = -Tr(p(0)F(t)A) (2.4) 

where 

F(t) = (1 -P)e'^*(l -P)L (2.5) 

This kernels is sufficient to determine the correlator of the noise in term of 
the n - point functions of the initial state since 

{FA{t)Fs{t')) = Tr {p{Omt)A){F{t)B)) (2.6) 
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The free field We have found more convenient to use the Fock represen- 
tation for the free field. Our basic objects are the creation and annihilation 
operators in real space a(x) and a^(x) satisfying 

[a(x),a^(x')] =(5(x-x') (2.7) 

They are related to the standard operators in momentum space by 

with dk = . The normal- ordered Hamiltonian then reads 

H^^J dxdx'd\:ii)h{^, x')a(x') (2.9) 

with 

/i(x,x') = j dke-'^^^'-'^'^cok (2.10) 
while the Heisenberg - picture operator 

a(x, t) = e'^'d{x)e-'^' = j dx'U{x - x'; i)a(x') (2.11) 

in terms of 

C/(x - x'; t)^ J dke-'^^^'-'^'^e-''"^' (2.12) 
2.2 The coarse - graining operator 

Let us assume a cubic lattice with side length L on a spatial hypersurface 
E of Minkowski spacetime. The centers of the cubes are identified with the 
points, coordinated by (nL, mL, rL) with n,m,r G Z. Let us denote, by 
[.] : S ^ S the function that takes a point x and assigns it at the point [x] of 
the center of the cube of the lattice in which x belongs. Also let /(x) denote 
the cube of the lattice in which x belongs. 

Our aim is to find a coarse graining operator, that intuitively does the 
following: to each field in spacetime it assigns a number of operators, which 
correspond to the smearing of the field over a lattice cube and assumed to 
lie on its centre . It is rather easy to construct this kind of operator by 
considering the smeared fields. 

Pa(/)=a(r) (2.13) 
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where is defined as follows 



r(x) = -^/ /(x')x/{x)(x')rfx' (2.14) 

where Xs(x) is the characteristic function of the set S. It is easy to check 
that this definition produces an indempotent, self adjoint with respect to the 
Hilbert - Schmidt inner product operator on the space of observables Q 
In a formal sense one can write the operator acting on unsmeared fields 

Pa(x) = ^ J dx'P{x, x')a(x') (2.15) 

where 

P(x,x') = x/{x)(x') (2.16) 
If one then employs equation (2.7) we obtain the interesting result 

Pa(x) = / (^^e"^'''''«(k)PL(k) (2.17) 

in terms of the function 

, sm(kiL/2) 

/..(k)=n4^ (2.18) 

This probably makes more transparent the meaning of the particular coarse - 
graining. The coarse grained field has a modifiacation in its spectral density 
and is essentially evaluated in the centers of the cubes in the lattice. It also 
implies the useful expression 

P(x,x') = J rffce-^''(M-'''VL(k) (2.19) 



^ This projector has a natural geometric significance which is worth stating. Recall 
that a field is an operator - valued distribution on Co(S), where Co(S) is the space of 
complex - valued functions of compact support on E. Now, the map [.] is a projection 
operator, projecting S to the lattice E/L'^. It is easy then to confirm that P is obtained 
from [.] by two successive applications of the pull - back operation. 
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3 The evolution equations 

Equipped with the above expression for the coarse graining operator, our task 
is now to write down an exphcit form of the equation (2.3). This involves 
the evaluation of the noise term and of the kernel 

G(r) = PL(1-P)e^^^(l-P)LP (3.1) 

We find it convenient to use the index notation employed in ||^, where an x 
dependence of a kernel is denoted by an index and integration by contraction 
of the indices (for more details see Appendix A of this reference). 
It is easy to verify that 

G(r)a'^ = H\U{T)\H',a'' (3.2) 

where 

H\ = P\h''dQ\ (3.3) 

Also Q = 1 — P, /i is defined by equation (2.9) and a bar denotes complex 
conjugation. Actually 

i/(x,x') = j dkdk'p{k)e-'^^''^+'^''''uJk 
X ({2nY6ik - k') - PL(k')A(k', k)) (3.4) 



where 



A(k, k') = J dye''^^y^*'^'y 



-Mkonl^ll^ (3.5) 

Substituting this together with (2.11) in the above equation we obtain for 
the kernel G 

G'(x,x';r) = / rfA;c/A;'dA;"e-^''M+*^"[^']wfcWfc/pL(k)pi(k')e-^^x'" 
X {{2nf6{k - k') - PL{k')A{k', k)) ((27r)35(k" - k') - pL(k')A(k', k"))(3.6) 
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The same steps would hold for the calculation of the noise terms. The 
corresponding kernel is easily found to be 

X ((27r)3<5(k - k') - pL(kO A(k', k)) {{2'n f6{\^ - k") - p^lk") A(k", k'));3.7) 

Now this equation together with its conjugate for a) are sufficient to deter- 
mine the hydrodynamic behaviour of the expectation values of the smeared 
fields fields. (Recall that free evolution does not mix the n-point function 
levels). Denote the coarse- grained Hamiltonian kernel by /i(x — x') . 

/i(x,x') = j dA;dA;'p(k)e-^^M+^'*'"'a;fePL(k')A(k',k) (3.8) 

In terms of the above kernels the relevant evolution equations are: 

i^a{pi.,i)— [ dx'h('K,x')a(x',t)+i [ drG(x, x'; r)a(x', t— r) = ^(x, t) (3.9) 
ot J Jo 

and its complex conjugate. ( We have denoted a(x) = Tr (p(O)Pa(x))). The 
noise can be modelled by the knowledge of its correlation functions. For 
example: 

(e(x,t))= y'rfxF(x,y)(a(y))o (3.10) 

(e(x,i)e(x',0) = / dydy'F{^,y;t)F{^',y';t'){d{y)d{y'))o (3.11) 

(e(x,t)e(x',t')) = / WF(x,y;t)F(x',y';t')(a(y)at(y'))o (3.12) 

(3.13) 

and continuing for the higher order correlation functions. 

4 The Gaussian approximation 

So far our approximations have been purely formal in the sense that no 
physical input has entered our equations. They have all been derived by 
foUwing the mathematics implicit in the Zwanzig technique. We still need to 
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enforce the three conditions we have stated in the introduction, in order to 
get meaningful hydrodynamic equations. 

The first observation is that the expectation value of ^ (x, t) is dependent 
solely on the expectation value of the field at the initial moment of time. 
Now, the scale by which (a(x)) varies is exactly what should be identified 
with the length l^ic- Assumption (2) implies that with good approximation 
we can set is to zero. Hence 



(e(t,x))^0 (4.1) 

Note that this condition does not imply, that the initial state has no irrelevant 
component. What is rather implied, is that the part of the state, where the 
one - point functions have support lies in the eigenspace of P. In particular 
higher order correlation lengths in the initial state can be larger than L and 
these are the ones that give rise to the noise 0. Note, that this assumption 
for the initial condition essentially fixes the origin of the time axis. 

Now, the difficulty of obtaining simple hydrodynamic equations is math- 
ematically connected to two facts: 

1. The complicated expression for the function p^. 

2. The appearance of the [.] in our exponential, which also gives rise to the 
presence of the function A. 

What we would like to have is to substitute [x] with x in the exponential and 
Pl with another easier to handle function. The condition for the former is 
essentially the possibility of taking the continuum limit for the lattice, which 
essentially amounts to the condition L « lobs- And we shall show that the 
later substitution possible by the use of the other hydrodynamic assumption 

The idea is instead of using the characteristic function for the set /(x) in 
equation (2.14) to use a smeared characteristic function. The most efficient 
way to perform this is to substitute the kernel P(x, x') with the Gaussian 

^(-'-')=(^^e-(— ')^/-^ (4.2) 



Since a quantum field is a system with an infinite number of degrees of freedom, in any 
state tliere is a priori an infinite number of lengtli parameters, subsets of them are roughly 
governing the behaviour of n-point functions. Even a Gaussian state is characterized by at 
least two such scales: the typical length scale of its mean value and the correlation length 
read from its two - point function, which are generically quite distinct. 
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Note that this substitution involves essentially two steps: 

1. Using a smeared characteristic function, which can be justified only if the 
width is much larger than the scale of variation for the expectation value of 
the field field, hence Imic « L. The error will be then of the order {Irnid LY- 

2. Substituting [x] with x, which is equivalent as we explained to the invo- 
cation of a continuum assumption. 

At this point we should stress that P(x, x') does not correspond to a true 
projector, since it is not indempotent. But in the particular subspace of the 
field states, that the hydrodynamic conditions of separations of the lengths 
is satisfied, it is very close to a true projector, with an error of the order 

m.iiy.{{lraic/ LY,L/lobs]■ 
Y{e^[lce equation (2.17) can now be writen in the form 

where 

p^k) = e-^'-^'/^ (4.4) 

and the function A becomes approximately 

A(k, k') = (27r)=^5(k - k') (4.5) 

Hence the kernels governing the evolution of the smeared fields simplify sig- 
nificantly 

G(x,x';r) = / dke-^'^^--'^ujl[{p~p'mfe-'-^^ (4.6) 

F(x, x'; t)^ j dke-'^^''-'''^u;ke-^'''p{k) (1 - p{k)f (4.7) 

^(x,x') = J dke-'^^''-'''^ p\k)ujk (4.8) 

Now the kernels G and h act on the coarse grained expectation values of the 
field, in which the high energy modes are suppresse by the Gaussian smearin. 
Hence within an approximation of {L/lohsY keep only the lowest order 

in the expansion of with respect to L. This gives 

^(x, x') = /i(x, x') (4.9) 

G'(x,x';r) = j fi^k^t^^g-ik(x-x')g-ia;fer 

= ^V^(m2 + V2)[/(x,x',r) (4.10) 



13 



On the other hand F acts on the unsmeared fields, hence we cannot drop 
the exponential damping term. But using the same argumentation we can 
expand (1 — pi) to obtain 



F(x,xV) = / 



^^g-ik(x-x'), 



ujke 



L2kV4g-L2k2 



V2t/(x,x') 



(4.11) 



with 



C/(x,xO= / 



^^g-ik(x-x'), 



cuke 



L2k2/4 



(4.12) 



which is essentially the relevant part of the evolution operator. Hence the 
kernels giving rise to dissipation and noise can be writen in terms of the full 
and the relevant part of the evolution operator. 

5 Remarks 

Our main results is essentially the Langevin type of equations (3.9) - (3.12) 
written in terms of the kernels (4.10) and (4.11). They should form the basis 
of construction of hydrodynamic equations, since having assumed classicality 
one can use the classical expressions to write energy and momentum density 
in terms of a(x) and a*(x). 

The equations admit simplifications in various regimes. For instance the 
dissipation term for a massless field reads simply 



where the non- local in time behaviour is encoded in the redarded time de- 
pendence of a. But, it has been impossible for this choice of coarse - graining 
to obtain a local in time equation, which could lead a field counterpart of 
the Navier - Stokes equation. 

One reason for this is probably to be found in our choice of coarse grain- 
ing. While physically realistic, it is conceivable that it might not be the most 
natural choice. The formalism of generalised quantum mechanics allows in 
principle space-time coarse grainings and it is probably through them that 
a maximal quasiclassical domain can be constructed. The possibility that a 




(5.1) 
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temporal coarse - graining is needed as well, is reinforced by examining stan- 
dard non- relativistic applications of the projector technique. The timescale 
of observation ought to be taken much larger than the characteristic time 
parameter appearing in the kernel G, if we want to obtain local in time (or 



Markovian) evolution equations |TT 



As far as the noise is concerned, it is highly unlikely that it can be approx- 
imated by white noise, even though for particular initial state the Gausian 
approximation might be a reasonable assumption. One needs to make phys- 
ical assumptions for the initial state, particular to the problem in hand 0. 

There are a re three interesting points to make, concerning our approach 

1. For our particular class of initial states, the contribution of noise and 
dissipation is very weak (the coupling constant is essentially L"^). Hence, 
while for short times one can safely rely on the evolution governed by the 
self dynamics of the coarse - grained field, as time increases predictability 
is going to be degraded, the faster the larger the long- range correlations in 
the initial state are. One would eventually have to use statistics (classical), 
constructing equations for the probabihty densities of the smeared field. 

2. The parameter governing both the loss of predictability and the flow of 
energy due to dissipation is the phenomenological coarse -graining parameter 
L. It is one of the great appeals of the histories approach, that they give an 
algorithm (computing maximal decohering histories), following which it can 
be determined. In our considerations we have encountered a feature, that 
could potentially augment our understanding of the quantum to classical 
transition: The parameters determining minimum coarse - graining such that 
decoherence is possible, will appear among the coupling constants determining 
the strength of dissipation and dijfusion in any hydrodynamic description 

3. In our coarse graining, essentially the low energy modes contribute to the 
relevant observables, but high energy modes (even if not directly coupled to 
the low energy modes) induce a significant noise in the equation of motion, 
by virtue of possible long - range correlations of the initial state. 

We should finally remark, that a more promising approach towards in- 



^ Note that in the projector technique an ignorance interpretation is assumed for the 
initial density matrix, that is the relevant component only can be known and the irrelevant 
should be modelled, usually by conjunction of physical arguments and application of the 
minimum entropy principle. While this is not a problem for most physical systems, in 
cosmology (where there is no meaning of an ensemble of systems) one has to be extremely 
careful on the assumptions involved for the modeling of noise 
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vestigating spacetime coarse grainings, implemented directly at the level of 
the field 's generating functional is currently under investigation. 
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